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O ■ I. INTRODUCTION 

o 

Historically, the spin-1/2 weakly interacting massive particles (WIMP) were considered as the first 
cold dark matter (DM) candidates. They interact with ordinary matter predominantly by means of axial 
vector (spin-dependent) and vector (spin-independent) couplings. 

Nowadays, the main effort in the direct dark matter search experiments is concentrated on the study 
of the spin-independent (or scalar) interaction of the dark matter particles with nuclei. It is due to a 
strong (proportional to the squared mass of the target nucleus) coherent enhancement of the dark matter 
particle scalar interaction with nuclei. The results obtained in the field are presented in the form of the 
^ ' exclusion curves for the total even rate as a function of the mass of the dark matter particles. The values 

of the cross section associated with the elastic scattering of WIMP due to scalar-nucleon interaction, 
which lie above these curves, are excluded. There is also the so-called DAMA contour which corresponds 
^J , to the first claim for evidence of the dark matter signal [lj . 

The main goal of this review is to attract attention back to the spin-dependent (or axial-vector) 
interaction of dark matter particles with nuclei. The importance of this type of interaction of the DM 
particles is due to the reasons as follows: i) the spin-dependent interaction of the DM particles provides 
us with twice stronger constraints on the SUSY parameter space in comparison with the spin-independent 
interaction, ii) in the case of spin-dependent interaction of heavy WIMPs with heavy target nuclei the 
so-called long g-tail behavior of the relevant form-factor allows detection of large nuclear recoil energy 
Oh, due to some nuclear structure effects; iii) it is worthwhile to note that by relying only upon the scalar 

interaction of the DM particles, which seems to be strongly suppressed, one might miss a DM signal Q. 
However, by a simultaneous study of both spin-dependent and spin-independent interactions of the DM 
particles with nuclei the chance for observing the DM signal is significantly increased. 

There are many different nuclear structure calculations (including the case of non-zero momentum 
transfer) for spin-dependent neutralino interaction with various nuclei, in particular with helium 3 He [2j, 



3 



O 



43 
Oh 



^ 



X 

5-H 



fluorine 19 F 0,1111, sodium 23]Nfa EHHHl. aluminium 27 A1 0, silicon 29 Si HHI1, 35C1 H.JEotassium 



nuorme "£ |^,U,DJ, sodium " u iMa Li |J. [J. H| . aluminium "Ai pj, silicon — oi jj,y,ia, ""<Ji |a|, potassium 
39 K 0, germanium 73 Ge [10, niobium 93 Nd H3, iodide 127 I @, xenon 129 Xe 0, 131 Xe BEHU, 
tellurium 123 Te [ll|, tellurium 125 Te [g, lead 208 Pb |3,|lj|. The zero-momentum transfer limits (mostly 
quenching) were investigated for the target nuclei Cd, Cs, Ba and La in [ill ll4L fl5J| . 

A dark matter event is an elastic scattering of a relic neutralino x (or X) ° n the target nucleus (A, Z), 
which results in a nuclear recoil with energy E^ detected by a proper detector. The differential event 
rate in respect to the nuclear recoil energy is a subject of experimental measurements. It depends on the 
distribution of the relic ncutralinos in the solar vicinity f(v) and the cross section of neutralino-nucleus 
clastic scattering |l(|, LLZl lla, Il9l 12(1 l2ll I22L |23| . The differential event rate per unit mass of the target 
material takes the form 

l!R =N^ f™dvf{v)v^{vrf). (1) 
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Here, N = N /A is the a number density of target nuclei. N and A stand for the Avogadro number and 
the atomic mass in AMU, respectively. The typical value of the nuclear recoil energy E-r, — q 2 /(2Ma) is 
about 10 _6 m x . Ma denotes the nuclear mass. 

The neutralino-nucleus elastic scattering cross section for spin- non-zero ( J ^ 0) nuclei is a sum of the 



coherent (spin- independent, or SI) and axial (spin-dependent, or SD) terms |8l Il3. I24J] : 
^ A . {v ^) = SI-MI 2 _ ^d(? 2 ) , S^) 



dq 2 ^^' nv 2 (2J+l) v 2 (2J+l) v 2 {2J+l) 

r sV0) F 2 D(g2) + ^ F|i(g2) . (2) 
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The normalized-to-unity (-Ffrj si(^) = 1) nonzero-momentum-transfer nuclear form- factors 

p2 r 2x _ gsr^si(g ) , 

^sd,siI<? J - ~Fa TTyT w 

°SD,Sn u J 

can be expressed through the nuclear structure functions as follows la, Il2l l24| : 

5 s A i(g) = E I<^I|Cl(<z)||J)| 2 ^|<J||Co(<z)||J)| 2 , 

L even 

SL(q) = E (KAr||T^ 5 ( g )||JV)| 2 + |(7V||£|(g)||7V)| 2 ). (4) 

Lodd 

Here, the double vertical lines denote the reduced matrix element. The explicit form of the transverse 
electric T el5 (q) and longitudinal C 5 (q) multipole projections of the axial vector current operator, scalar 
function Cl(<?) and Sgj grj(9) at zero momentum transfer can be found in Appendix A. For q — the 
nuclear SD and SI cross sections (in (J2J) take the forms 



4/4Ssi(0) _/4 „2_p 
(2J+1) "M 2 SI 



^i(O) = ToV, iV =Tf^Xi(0). (5) 



^n(O) = ^jf ^^K(S-) + an (S-)} 2 (6) 

= ^r-^T { V^d(OXS^) +sign(a p a n ) v /^ D (0)(S^)J . (7) 

Here, ^ = is the reduced mass of the neutralino and the nucleus and it is assumed that 

m x + Ma 

H„ = jj/z. The dependence on effective neutralino-quark couplings C q and A q in the underlying (SUSY) 

theory (see Appendix A) 

C eff = E (A ' X7m75X • <77 M 75<7 + C q ■ XX • 9<z) + ■■• (8) 

and on the spin (A q p ' n ') and the mass (fq ) structure of nucleons enter into these formulas via the 
zero-momentum-transfer proton and neutron SI and SD cross sections: 

4M = 4^4 cr = E c ^ (p ' n) ; (9) 

7T * — ' 



o*£(0) = 12^X„ « P = ^^ p) , «-E^ n) - (10) 



The factors A q p ' , which parameterize the quark spin content of the nucleon, are defined as 2Ay 1,p )s M = 
(p, s|^ g 7 A '75V'(j|p, s)( P ,n)- The total nuclear spin (proton, neutron) operator is given by 

A 

i 

where the index i runs over all nucleons. 

The expectation values of the spin and angular operators are evaluated, as a rule, in their z-projection 
by assuming the state with the maximal value of the angular momentum projection Mj = J: 

(S) = (N\S\N) = (J, Mj = J\S Z \J, Mj = J). (12) 



Thus (S p („)) is the total spin of protons (neutrons) averaged over all nucleons of the nucleus (^4, Z). 

The mean velocity (v) of the relic neutralinos of our Galaxy is about 300 km/s = 10~ 3 c. Assuming 
<7max-R *C 1, where R is the nuclear radius and q m ax = ^f-AV is the maximum of the momentum transfer in 
the process of the yA scattering, the spin-dependent matrix element takes a simple form (zero momentum 
transfer limit) |a,l3 : 

M = C(N\a p S p + a n S n \N) ■ s x = CA(N\J\N) ■ s x . (13) 

Here, s x denotes the spin of the neutralino, and 

(JV|q p S p + o w S n |JV) (JV|(fl p S p + fl n S n )-J|iV) 

(N\3\N) J(J+l) ' ( ' 

Note a coupling of the spin of \ to the spin carried by the protons and the neutrons. The uncertainties 
arising from the electroweak and QCD scale physics are incorporated in the factors a p and a n . The 
normalization factor C involves the coupling constants, the masses of the exchanged bosons and the 
mixing parameters relevant to the lightest supersymmetric particle (LSP), i.e., it is not related to the 
associated nuclear matrix elements J2"| . The above conclusions concerning the spin-dependent part of the 
neutralino-nucleus scattering amplitude are also valid for the amplitude of any Majorana WIMP-nucleus 
scattering process. In the limit of zero momentum transfer q = the spin structure function in Eq. J3J 
reduces to the form 

Siu(0) = — L - A 2 J(J+1). (15) 

The nuclear matrix element A4 in Eq. I|13l) is often related to the matrix element of the nuclear magnetic 
moment, which also consists from the matrix elements of the total proton and neutron spin operators: 

H = (A%*S n + g l n L n + g ;S p + g l p L p \N). (16) 

The free particle g-factors (gyromagnetic ratios) are given (in nuclear magnetons) by: 

£ = -3.826, g l n = 0, g ; = 5.586, g l p = 1. (17) 

The nuclear magnetic moment /i is often used as a benchmark for the accuracy of the calculation of S p 
andS„ [EH- 

If the neutralino mass m x is larger than few tens of GeV, the value of the product qR is no longer 
negligible and the so-called finite momentum transfer limit has to be considered in the case of the neu- 
tralino scattering on medium-heavy and heavy nuclei. The corresponding formalism is a generalization of 
that used for the description of weak and electromagnetic semi-leptonic interactions in nuclei. We shall 
follow the conventions of @, 13 • There is an advantage to use the isospin instead of the proton- neutron 
representation when discussing %A scattering at finite momentum transfer. By rewriting the isoscalar 
and isovector coupling constants as a§ — a n + a p and a± = a p — a n , respectively, the spin-dependent cross 
section S^ D (q) decouples into the isoscalar term Sqq, the isovector term Sn and the interference term 
Soi as follows: 

Ssd(«) = aj}Sbo(g) + a\ S n (q) + a aiS i(q). (18) 

The structure functions in Eq. I|18(l consist of the expectation values of operators JL{o_r)\Yi J o-\ L±1 (L 
even), which depend on spin and spatial coordinates. Using the decomposition of Sa (q) in l|18|) one can 
obtain structure functions for a x of arbitrary composition. 

The cross section of neutralino-nucleus scattering at zero momentum transfer exhibits a strong depen- 
dence on the details of the nuclear ground state p. The goal of this review is to collect the results of 
different calculations and discuss their spread and relevance. 

II. THE ZERO MOMENTUM TRANSFER LIMIT 

Only nuclei with odd number of either protons or neutrons possess nonzero total nuclear spin. At 
first, the independent single-particle shell model (ISPSM) was employed by Goodman and Witten |26j | 
and later by others [H, 123 u3 to estimate the spin content of the nucleus for the detection of dark 
matter. This model utilizes the shell structure of the nucleus, in particular the fact that if certain magic 
numbers of nucleons occur in the nucleus, it exhibits remarkable stability properties, and the ground 
state expectation values of the total spin J and the parity of the nucleus can be described by those of 
the extra nucleoli. 



For nuclei whose angular momentum J is given by a single neutron (proton) with spin s and the orbital 
momentum L, J = L + s, (the even number of nucleons which remain form pairs with the opposite angular 
momentum projection and zero spin, i.e., they do not contribute to J, we have |28j : 



_ J(J+l)-L(L + l) + l 

n(p)/ 2 J + 2 ' P{n)l 



K M ) = - hrr^ — > (Km)=V- (19) 



In the ISPSM the entire angular momentum J and the parity of the odd-odd nucleus (A,Z) are identified 
with a single proton (Z-odd) or neutron (Z-even) state. Then the spin matrix elements are given by Eq. 
P)l. 

The ISPSM offers only a rough estimate of the spin matrix elements. From Tables ITl lXIIII it follows 
that the ISPSM predictions significantly overestimate the values obtained in realistic calculations. The 
ISPSM results are qualitatively good only for light nuclei with the single nucleon outside the closed shell 
(e.g., 17 0); however, they become increasingly poor for heavier isotopes, especially for those with many 
particles outside the closed shells. The realistic calculations take into account the complex structure of 
the nuclear wave functions, the fact that the contributions to spin matrix elements both from paired 
nucleons and unpaired ones cannot be neglected and the phenomenon that the free nucleon structure 
coefficients are renormalized when nuclear medium effects are relevant. These effects arc known to play 
an important role in the calculation of the matrix elements of the magnetic dipole moment too (see, for 
example, |15J]). 

There are several elaborated nuclear structure approaches which lead to more accurate predictions 
of spin matrix elements associated with the dark matter detection on nuclei in comparison with the 
ISPSM. A full list of these models,to our knowledge, includes the Odd Group Model (OGM) J2^ 
and the extended OGM (EOGM) J2J, l2p o f Engel and Vogel, the Interacting Boson Fermion Model 
(IBFM) of Iachello, Krauss, and Maino [l5J], the theory of Finite Fermi Systems (TFFS) of Nikolaev 
and Klapdor-Kleingrothaus [2£J, the Quasi Tamm-Dancoff Approximation (QTDA) of Engel [l2 |. the 
nuclear shell model (SM) applied by Pacheco and Strottman [lj|, Engel, Pittel, Ormand and Vogel [Hi ]. 
Engel, Ressell. Towner and Ormand p|, Ressell et al. [§], Ressell and Dean [fg; and by Kosmas, Vergados 
ct al. 0,0, uM t° different nuclear systems, the so-called "hybrid" model of Dimitrov, Engel and Pittel 
[jj and the perturbation theory (PT) based on calculations of Engel et al. p| • 

The ISPSM predictions are fairly accurate for near-closed-shell-nuclei, but further away they tend to 
overestimate the spin contribution to the magnetic moment. In an open-shell nucleus, the last odd particle 
polarizes the other nucleons in the direction opposite to its own spin, which results in a spin-quenching 
effect entirely absent in the single-particle picture. Denying the idea about importance of only the last 
odd nucleon Engel and Vogel arrived at the "odd-group" model 29] by assuming that the nuclear spin is 
carried by the "odd" unpaired group of protons or neutrons and only one of either (S„ ) or (Sp ) is non- 
zero. The odd-group spin matrix elements are expressed with the measured nuclear magnetic moment fi 
as 



V - 9p J »-J ,cA\ i l - 9 l n J 
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r> gs-gl p 4.586' { * n} 5 * - 5 < 3.826' l Uj 



where g denote gyromagnetic factors of the free nucleon [see IJ17JI ]. The OGM has been found successful, 
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e.g., in the case of 29 Si with J = ~ and unpaired neutrons. The experimental value of the nuclear 



magnetic moment /i = —0.555 implies (S 29 ) ss 0, and (S 29 ) w 0.15, which is |lfj in good agreement 
with the shell- model calculation of Ressell et al. |j|. The results of OGM calculations are collected in 
Tables HV lXIIII We note that for 73 Ge with a complex nuclear structure the odd-group model prediction 
disagrees with the realistic calculation of 0, Oj ■ 

The odd-group model is a significant improvement in comparison with the ISPSM. The weak points 
of this approach are the facts that the roles of small but not vanishing angular momenta of the even 
system and of the meson-exchange currents, which can renormalize the g factors in Eq. H20|) . are ignored. 
Engel and Vogel improved the OGM [23] by using additional information about /3-decay ft values and 
measured magnetic momenta of "mirror pairs" for nuclear systems with (A < 50). For these nuclei they 
proposed to use two relations |29l ]: 

(^) 2 «s odd )-<s cvcn )) 2 = f^-i^ J 



gv \ ft J J+l 

Mis = J + 0.76((S odd ) + (S cvon »+^. (21) 

Here, gA and gv are the axial vector and vector coupling constants, respectively, n x is a small correction 
induced by heavy meson exchange and /i IS is a sum of two mirror magnetic moments (isoscalar moment). 
For free nucleons we have gA = 1-25 and gv = 1-0. However, in the nuclear matter due to the effect 
of renormalization the value gA/ 'gv = 1-00 ± 0.02 is often considered. For light nuclei the spin matrix 



elements (Sf) and (S^) evaluated within the extended odd group model (EOGM) [2!j are listed in 
Tables fll liVI We see that there is quite good agreement between the EOGM (with gA/gv = 1-00) and 
the more sophisticated shell-model calculations for light odd-even isotopes performed by Pacheco and 
Strottman 14]. Their calculations for A < 16 nuclei assumed the Cohen-Kurath interaction J31J and a 
complete basis within the p-shell model space. For A > 16 the Rcid interaction was considered and the 
basis consisted of all allowed states within the ls-Od shell-model space. The results of [l4j are given in 
Tables ITl lIVI For heavier mirror nuclei with A close to 50 the shell-model calculations are difficult due to 
a large amount of configurations which have to be taken into account. 

The light nucleus 27 A1 is one of the active ingredients of a very high-resolution and low-threshold 
sapphire-crystal (A^O^-based detector for dark matter search. Engel, Ressell, Towner and Ormand Q 
performed calculation of proton and neutron spin expectation values for this isotope with the help of 
the Lanczos m-scheme shell- model code CRUNCHER |32j . The nucleus 27 A1 lies in the middle of the sd 
shell and the m-scheme basis for 27 A1 contains 80115 Slater determinants. Good agreement between the 
calculated and measured spectroscopy of excited states was achieved for this nucleus. In addition, the 
experimental value of the magnetic moment /i cxp = 3.6415 /ijv (mjv = eh/2m p is the nuclear magneton) 
was reproduced in calculation well. The theoretical value /i = 3.584 fix was obtained with the help of 
the free-particle g- factors Q. We recall that calculation of the magnetic moment requires evaluation 
of the same spin matrix elements needed to determine the WIMP structure functions at q 2 — 0. The 
corresponding values of (S 27 ) and (S 27 ) are given in Table IIIII We note that the authors calculated 
structure functions S(q) of 27 A1 at q ^ as well Q. 

In the case of 39 K the shell-model diagonalization needed for the calculation of the nuclear spin matrix 
elements requires severe truncations to the active model space. The problem is that 39 K is so near the 
boundary between the sd and pf shells and excitations of particles into higher shells can have significant 
effects that are often not well simulated by effective operators. Thus, for this nucleus Engel, Ressell, 
Towner and Ormand |jj used an alternative scheme based on perturbation theory (PT) for the evalua- 
tion of spin matrix elements. It was successfully implemented in calculations of several spin-dependent 
observables in closed-shell- plus (or minus)-one nuclei |22| ■ The details of the method and the calculations 
can be fobbed in Ref. [3. The authors considered two different residual interactions. One (denoted as / 
in Table llVJl is related to the one-boson-exchange potential of the Bonn type, but it is limited only to four 
or five important meson exchanges. The resulting interaction has a weak tensor-force component typical 
of Bonn potentials. The other (denoted as 77 in Table HV|) is represented by full G-matrix elements of the 
Paris potential parameterized in terms of sums over Yukawa functions of various ranges and strengths. 
Interaction II exhibits a strong tensor force. The quality of the wave functions obtained was judged 
in terms of magnetic moments and Gamow- Teller matrix elements, including meson-exchange currents, 
isobar currents, and other relativistic effects. The results were presented for both isoscalar and isovector 
magnetic moments and their sum, i.e., for the magnetic moment of 39 K, whose value is given in Table llVl 
The magnetic moments calculated with the help of both interactions differ only slightly from each other 
and showed good agreement with corresponding experimental values. The same nuclear wave functions 
of 39 K were also used for the calculation of (S 39 ) and (S 39 ) ( Table llVJl and the structure function S(q) 
7] . It is worthwhile to notice that for both types of interactions the zero-momentum transfer spin matrix 
elements coincide well with those obtained within the phenomenological EOGM 29]. 

For dark matter targets constructed of heavy nuclei, in particular, Ge, I, and Xe, the first elaborated 
calculation of spin-dependent matrix elements relevant to WIMP scattering was performed by Iachello, 
Krauss, and Maino within the Interacting Boson Fermion Model (IBFM) [13. The applied IBFM 
wave functions were tested in a comprehensive analysis of excitation energies, electromagnetic transition 
rates and intensities of transfer reactions 15]. In this model the total spin operator has the form S = 
S s p + S s n + s pm where sT are the paired proton (neutron) spins, and s ,-. is the remaining 
unpaired proton (neutron) spin. To estimate the matrix elements of the paired nucleons one should 
know the structure of Cooper pairs (bosons), which were incorporated in the model by fitting the matrix 
elements of the magnetic moments. The authors end with the conclusion that the ISPSM predictions are 
generally within 15% of their results for the nucleon spin matrix elements 15] . If quenching of free-particle 
coefficients in the nuclear environment is considered, the value of the spin matrix elements is reduced 
by an additional factor not exceeding 60%. The results obtained are listed in Tables IVlllXIl We note 
that the spin-dependent WIMP-nucleus cross section associated with the IBFM spin matrix elements is 
always smaller than the ISPSM prediction but not more than a factor of 5 [jj|. While the IBFM can 
incorporate the dominant collective effects, it has some difficulty in including the spin polarization, which 
plays a crucial role in axial vector scattering. Unfortunately, this approach cannot be readily applied to 
the case of nonzero momentum transfer [5j . 

In J5J Vergados with co-authors investigated the spin-dependent elastic neutralino scattering with light 
nuclei 19 F, 23 Na, and 29 Si. The spin contribution to the differential cross section was obtained by 
the shell model calculations in the sd shell using the Wildenthal interaction, which was developed and 
tested over many years. This interaction is known to reproduce accurately many nuclear observables 



for sd shell nuclei. The Wildenthal two-body matrix elements as well as the single-particle energies are 
determined by fits to experimental data in nuclei from A = 17 to A = 39. The shell-model wave functions 
used by the authors were tested in the calculation of the low-energy spectra and ground state magnetic 
moment. Rather good agreement between theoretical results and experimental data was achieved. This 
fact increases the confidence level of the calculated spin matrix elements which are listed in Tables HTl II VI 
It is worth mentioning that these spin matrix elements are in good agreement with those of previous 
calculations [a| . The authors of Ref. Q found 19 F to be the most favorable target for dark matter search 
via spin-dependent interaction of relatively light dark matter particles. It is favored due to the fact that 
the corresponding spin matrix element is not quenched and that various isospin channels add coherently. 
Further, it was demonstrated that the effect of the nuclear structure on the elastic scattering cross section 
of LSP with light nuclei (including q ^ behavior) is well understood, both for coherent and spin modes 

Si _ 

The Theory of Finite Fermi Systems (TFFS) was used by Nikolaev and Klapdor-Kleingrothaus |30| 
to describe spin matrix elements in the nuclear medium and to evaluate quenching of zero-momentum 
nuclear spin matrix elements of heavy nuclei due to residual interactions. Contrary to the OGM and the 
IFBM studies the TFFS calculation of spin matrix elements is not related with the experimental value 
of the associated magnetic moment. The TFFS proton and neutron spin averages (S p / n \) are suppressed 
in comparison with the corresponding ISPSM predictions and in some cases they differ significantly from 
the OGM values too J2J|. However, they agree well (with exception of the case of 73 Ge) with the results 
obtained by Pacheco and Strottman 14] in a completely different IFBM approach fTablcs lVll 1X1(1 . 

The momentum transfer dependence of the structure function S(q) associated with scattaring of dark 
matter particles from 131 Xe, a promising heavy target for the dark matter search experiment, was in- 
vestigated by Engel [IJ by using the configuration-mixing quasiparticle Tamm-Dancoff approximation 
(QTDA). In the zeroth order the ground state of 131 Xe was represented as the ld 3 / 2 quasineutron excita- 
tion of the even-even core |0) treated in the BCS approximation (BCS-based model of the Fermi surface). 
In the case of odd-multipole operators (@J the one-quasiparticle approximation corresponds to the ISPSM 
approach of |28j . In order to incorporate nuclear structure corrections originating from the residual in- 
teraction three-quasiparticle configurations of the form [t / ^ 3 / 2 [^fe I/ /] K ] 3 ^ 2 |0) and [^ 3 / 2 [ 7r fe 7r j :]" R "] 3 ^ 2 |0) were 

admixed. Here iv and v 1 represent the proton and neutron quasiparticle creation operators, K is an 
arbitrary intermediate angular momentum, and k, I run over a valence space consisting of the 2s, Id, Og 
and Oh harmonic oscillator levels [l2| . Despite the fact that the amplitudes associated with the admixed 
three-quasiparticle states are small (less than 5%), these admixtures can lead to a substantial effect. The 
experimental value of the magnetic moment of 131 Xe, which is about 0.69 /l/jv, was reproduced with an ac- 
curacy of 2% in the QTDA (note that the ISPSM value is almost twice larger). The same approximation 
scheme results in (S^ 1 ) = -0.041 and (S r \ 31 ) = -0.236 (Table EJ. 

The neutralino-nucleus cross section associated with the spin-dependent interaction is determined by 
the distribution of the spin in the nucleus. This observable is difficult to describe accurately as the lowest 
order nucleons create pairs with zero spin. The spin-dependent scattering mostly takes place near the 
nuclear Fermi surface and it is affected by the behavior of relatively few nucleons. The silicon nucleus 
is a relatively light. Thus, shell-model calculations or arguments based on existing data of magnetic 
moments (OGM and EOGM) allow a reliable prediction of deviations from the simple ISPSM picture. 
A different type isotope is germanium with a complex nuclear structure. For this nucleus the reliable 
calculation of spin expectation values is rather difficult. Niobium isotope 93 Nb is something in between 
the above two special cases. It is a heavy nucleus, which can be represented by a basic shell-model 
space corresponding to three protons in the lpi/2 or Ogg/2 levels and two neutrons in the 1^5/2 level 
|10| . This model space was considered by Engel, Pittel, Ormand and Vogel [l(| in the calculation of the 
magnetic moment for this isotope. They found \i to be equal to 6.36 /l/jv, which exceeds the experimental 
value // = 6.17 [i at. In order to obtain better agreement with the experiment the authors considered a 
"large" model space which included all basis states in which one proton or one neutron is excited from 
the small model space. Then, they ended up with the magnetic moment equal to 5.88 /xjy, a value that 
is smaller. The explanation could be that meson-exchange currents renormalize orbital proton g-factors 
upwards by about 10%, increasing the /i without altering the values of the nuclear spin matrix elements. 
The discrepancy between the (Sp 3 „) results obtained within the small and large model spaces ( Table Ivij) 
provides an indication of uncertainty of the nuclear-structure calculation [l(| . 

Germanium isotopes (especially large-spin 73 Ge) are considered to be the most promising material for 
the direct dark matter search experiment. However, there are fundamental difficulties in describing, e.g., 
the spin content of 73 Ge due to its complicated collective structure. Several studies were devoted to nuclear 
structure aspects of spin-dependent scattering of neutralinos from 73 Ge. Engel and Vogel (OGM) [23 
used measured magnetic moments to estimate the quenching of the nucleon spin in several heavy nuclei, 
including germanium. Iachello, Krauss and Maino [jj| employed the IBFM, and Nikolaev and Klapdor- 
Kleingrothaus [3(j used the TFFS to calculate the same quantities. There are two most comprehensive 
spin structure analyses for 73 Ge. A large-basis shell model study was performed by Ressell et al. 0, 



who calculated the full spin-dependent neutralino response including g-dependence of form factors. An 
equally comprehensive calculation was realized by Dimitrov, Engel and Pittel 9]. The authors obtained 
significantly different results in comparison with other studies and they argue that their results are more 
reliable than the previous ones. 

By using a reasonable two-body interaction Hamiltonian and appropriate large model spaces Ressell et 
al. [3 calculated the ground state wave functions for 29 Si and 73 Ge. In particular, the universal sd shell 
interaction of Wildenthal was used for calculation of the wave functions for silicon. The nuclear wave 
functions obtained were tested in the analysis of energy pattern of excited states and magnetic moments. 
Once reasonable agreement among theoretical results and experimental data was obtained, the ground- 
state wave functions were used to calculate the neutralino-nucleus nuclear matrix elements. In addition, 
finite momentum transfer matrix elements and cross sections for the spin-dependent elastic scattering 
of neutralinos from 29 Si and 73 Ge Q were evaluated. The computations were performed by using the 
Lanczos method m-scheme nuclear shell model (code CRUNCHER) [33] ■ The m-scheme basis for 29 Si 
had a dimension of 80115 Slater determinants. In the limit of zero momentum transfer the scattering 
matrix elements for 29 Si are in general agreement with previous estimates (Table Hvjl . 

For the study of 73 Ge Ressell et al. [8j chose the Petrovich-McManus-Madsen- Atkinson interaction [34J , 
which is a reasonable approximation to a full G-matrix calculation. This interaction proved to be both 
adequate and tractable in shell model applications. Two different model spaces were considered. The 
"small" space was determined by an m-scheme basis dimension of 24731 Slater determinants. The "large" 
space allowed much more excitations with an m-scheme basis dimension of 117137 Slater determinants. 
Despite fairly large size of the bases, rather severe truncations in the space were enacted. The small 
space is the smallest one in which it is possible to obtain agreement with the experimental spectrum 
energy levels. The dimension of the large basis was limited by the computer time and the memory 
storing constraints [8j . No phenomenological interaction has been developed for Ge-like nuclei and fairly 
severe truncations to the model space have to be imposed to obtain manageable dimensions. Despite 
these obstacles, the ground state wave function for 73 Ge allowed good description of the low-lying excited 
states and the ground state to ground state spectroscopic factor. The large model space wave function 
of Ge led to an improved description of the ground state expectation values, in particular of the value 
of the magnetic moment, in comparison with the ISPSM and IBFM estimates. The calculated magnetic 
moment \i from Q exceeds the experiment value, but the authors stressed that the same quenching of 
both (i and the Gammov- Teller (GT) spin matrix elements was almost universally required in shell model 
calculations of all heavy nuclei. Assuming the isovector spin quenching factor to be 0.833, agreement 
with the measured fi is obtained. In principle, it is not obvious that quenching is really needed in 
neutralino- 71 Ge scattering but if so, Ressell et al. believed that the correct answer might be in the range 
between the quenched and unquenched values. It was found (Table ^TJ that the zero- momentum-transfer 
spin-neutron matrix element (S 73 ) of 73 Ge was a factor of 2 larger than the previous predictions (except, 
obviously, the ISPSM value). Thus, even if quenching is assumed, the calculated scattering rate is about 
twice as large as any of the estimates made before ii] ■ 

A different sophisticated approach for evaluation of the spin structure of Ge was considered by 
Dimitrov, Engel and Pittel. It relies on the idea of mixing variationally determined Slater determinants, 
in which symmetries are broken but restored either before or after variation. This approach is described 
in detail in [35j • In the calculation of [|| the symmetries broken in the intrinsic states are those associated 
with rotational invariance, parity, and axial shape. The hybrid procedure used restores axial symmetry, 
parity invariance, and approximate rotational invariance prior to the variation of each intrinsic state. 
Subsequently, before mixing the intrinsic states the rotational invariance is fully restored. The procedure 
allows fully triaxial Slater determinants at the expense of particle-number breaking. The results of [35j 
indicate that the trading of number nonconservation for triaxiality is a good idea, despite the apparent loss 
of pairing correlations traditionally associated with the former. Pairing forces evidently induce effective 
triaxiality. The numerical results J23 show that the approach is accurate and efficient for describing even- 
even systems while also providing reliable reproduction of the collective dynamics of odd-mass systems 

a 

For 73 Ge the calculations were performed by assuming, both for protons and neutrons, a single-particle 
model consisting of the full 0/, lp shell and the 0<?g/2 and 0<?7/2 levels. The main idea was to include 
all of the single-particle orbits that could play an important role in reproducing low-energy properties of 
the 73 Ge |9j. It is well-known that a crucial ingredient in any realistic nuclear-structure calculation is 
the appropriate form of the nuclear Hamiltonian. The one- and two-body parts of the Hamiltonian have 
to be compatible with each other as well as with the model space. This is difficult to achieve because 
microscopic two-body interactions, derived for example from a G-matrix, include monopole pieces that 
are unable to describe the movement of spherical single-particle levels as one passes from the beginning 
to the end of a shell. A proposed procedure for avoiding this problem consists basically in removing 
all monopole components from the two-body interaction and shifting their effects to the single-particle 
energies. This procedure was used by Dimitrov, Engel and Pittel [3 — their two-body force was a fit to 
Paris-potential G-matrix modified as just described above. The calculated ground-state magnetic dipole 



moment is in good agreement with the experimental value. Ressell et al. [8j in their large space shell 
model calculation were able to reduce /i significantly to — 1.24/1 ^ (without direct quenching) but could 
not account for the remaining difference. On the contrary, the calculation of Dimitrov, Engel and Pittcl, 
despite the small number of intrinsic states, contains the full quenching required by experiment [9|. By 
making a comparison with the results of Ressell et al. again || , significant disagreement is found for the 
neutron spin. The calculated value of Dimitrov, Engel and Pittel is significantly smaller (Table Ej . The 
large and negative neutron spin g-factor (g^ — —3.826) is favored by the correct \x value. The differences 
in the spins, unlike those in the orbital angular momenta, carry over into WIMP scattering cross sections. 
Thus, following Ressell et al. Q, no significant increase is expected in the neutralino- 73 Ge scattering rate. 

The advantage of Dimitrov, Engel and Pittel's approach for calculation of neutralino cross sections 
is that it correctly represents the spin structure, requires neither quenching at q — nor arbitrary 
assumptions about the form factor behavior at q ^ Q- The spin matrix elements depend in general 
rather sensitively on the details of the nuclear structure. Since the matrix elements at q = are often 
quenched, the momentum dependence of the matrix elements was more important than it was naively 
expected. As a matter of fact, one has to include a lot of configurations to accommodate all multipoles, 
which result in very large Hilbert spaces in complex nuclei like 29 Si and 73 Ge. It will be therefore a very 
hard task to substantially improve the calculations of Ressell et al. [g and Dimitrov, Engel and Pittel 
[jj for these elements. 

For evaluations of the spin matrix elements in the heaviest possible nuclei relevant to dark matter 
search Kosmas and Vergados have chosen 207 Pb 0, Il3- Among the targets which were considered for 
direct neutralino detection, 207 Pb stands out as an important candidate. The spin matrix element of this 
nucleus has not been evaluated quite accurately, since one expected that the neutralino spin interaction 
is important only with light nuclei. But the spin matrix element in the light systems is quenched. On 
the other hand, the spin matrix element of 207 Pb, especially the isoscalar one, does not suffer unusually 
large quenching, as is known from the study of the magnetic moment. It is believed that Pb has 
a quite simple structure, its ground state can be described as a 2pi/ 2 neutron hole outside the doubly 
magic (closed-shell) nucleus 208 Pb. Due to its low angular momentum, only two multipoles L = and 
L = 2 can contribute even at large momentum transfers. One can thus view the information obtained 
from this simple nucleus as complementary to that of 73 Ge, which has a very complex nuclear structure 
0, Hil- -"- n the q = limit Vergados and Kosmas gave the spin matrix element in the simple form 

|J| 2 = /°fi (0) + /|Oi(0) , and found that Q o (0) = -0.95659/V3, and ^i(O) = 0.83296/^3 0,113. 

These values were recalculated in the form of spin variables (S p ( n )) given in Table IxTTTl 

Ressell and Dean have performed most accurate nuclear shell model calculations of the neutralino- 
nucleus spin-dependent or axial cross section for several nuclei in the A = 127 region, which are important 
for dark matter search. Their set of structure functions S(q) is valid for all relevant values of the 
momentum transfer. Conventional nuclear shell model of Wildenthal 36] quite accurately represents 
spin-dependent neutralino-nucleus matrix elements when a reasonable nuclear Hamiltonian is used in a 
sufficiently large model space [(j. Until recently, both of these ingredients have been absent for nuclei 
in the ?>s2dlg-j /2I/111/2 shell. Ressell and Dean considered two residual nuclear interactions based upon 
recently developed realistic nucleon-nucleon Bonn A 37) and Nijmegen II [38j potentials. These two 
nucleon-nucleon potentials were used in order to investigate the sensitivity of the results to the particular 
nuclear Hamiltonian. 

The Bonn-A-based Hamiltonian has been derived for the model space consisting of the 
lfl'7/2) 2^5/2, 3s]/ 2 , 2d 3 / 2 , and l/in/2 orbitals, allowing one to include all relevant correlations. In or- 
der to get good agreement with observables for nuclei with A ss 130, the single-particle energies (SPEs) 
were adjusted. The SPEs were varied until reasonable agreement between calculation and experiment 
was found for the magnetic moment, the low- lying excited state energy spectrum, and the quadrupole 
moment of 127 I. Once the SPEs are specified, a reasonable Hamiltonian can be used for the nuclei under 
investigation. 

To perform a full basis calculation of the 127 I ground state properties in the space consisting of the 
157/2, 2<i5/2, 3si/2, 2<i3/2, and lft.11/2 orbitals, one would need basis states consisting of roughly 1.3 x 10 9 
Slater Determinants (SDs). Current calculations can diagonalize matrices with basis dimensions in the 
range 1-2 xlO 7 SDs. Therefore clearly severe truncations of the model space are needed J6J. Fortunately, 
given the size of the model spaces that can be treated, a truncation scheme that includes the majority 
of relevant configurations can be devised. Finally (after relevant truncations, see details in [y]) the m- 
scheme dimension of the 127 I model space is about 3 million SDs. The calculated observables agree well 
with experiment. These interactions do not seem to prefer excitation of more than one extra neutron 
pair to the l/in/2- Most configurations have six neutrons in that orbital, while eight are allowed. Hence, 
this model space is more than adequate. It is this truncation scheme that was used for the two Xenon 
isotopes considered (A — 129 and 131). 

In almost every instance, the results of 6] ( Tables IVlIIllXll show that the spin (S^) (i = p, n) carried 
by the unpaired nucleon is greater than that found in the other nuclear models (except for the ISPSM, 



where (Si) is maximal). Despite these larger values for (Si), these calculations have significant quenching 
of the magnetic moment and are in good agreement with experiment in all cases. The larger values 
of (Si) are due to the fact that more excitations of the even group of the nuclei were allowed p|. The 
differences in the response due to the two forces is clearly visible in Tables IVIIII Ixl In all cases reasonable 
agreement between calculation and experiment for the magnetic moment (using free particle (/-factors) 
is achieved. It is obvious that the differences between the two calculations are non-trivial but they are 
quite a bit smaller than the differences coming from the use of alternate nuclear models. This shows that 
the interaction is not the primary uncertainty in calculations of the neutralino-nucleus spin cross sections 

@ 

The results obtained by Ressell and Dean give a factor of 20 increase in iodine's sensitivity to spin- 
dependent scattering over that previously assumed. Due to the form factor suppression a sodium iodide 
detector's spin response is still dominated by 23 Na but not to the extent previously thought. For the 
remainder of the nuclei considered Tables lVIIIhlXl also reveal increased scattering sensitivity, though much 
more modest gg. 

Before finishing this section we, following Ressell and Dean 6], discuss the quenching problem and some 
related uncertainties. As is already noted above, the comparison of the computed magnetic moment and 
its experimental value has been used as the primary indicator of the calculation's reliability. This seems 
quite reasonable in the light of the similarities between the matrix elements in Eqs. (|13|l and l|16(l . This 
prescription is not free of several potential problems 0, Q . Not only does fj, depend upon the orbital 
angular momentum Li but the spin angular momentum Si is subtly different. The neutralino-nucleus 
matrix element Ijl3(l results from the non-relativistic reduction of the axial- vector current. Because of this, 
it is not strongly affected by meson exchange currents (MECs). The magnetic moment's spin operators, 
Si, are a result of the non-relativistic reduction of the vector current. They can be strongly affected by 
MECs 7] . The effects of MECs upon /i is typically lumped together with several other effects to give 
effective g-factors. Unfortunately, there is no hard and fast rule as to what effective ^-factors are the best. 
One usually chooses to remain with the free particle (/-factors. As an example of the potential uncertainties 
this ambiguity leads to, the calculated magnetic moments for these nuclei based on a reasonable set of 
effective (/-factors were also included in Tables IVIIIllXl The "quenched" magnetic moments are the 
values in curled parentheses and the effective g-factors used are g^ = —2.87, g l n = —0.1, gt, = 4.18, and 
g l p = 1.1. The tables show that these (/-factors do little to improve the concordance between calculation 
and experiment 6]. A related concern involves the quenching of the (isovector) Gamow- Teller (GT) 
g-factor, gA 0, la|- The spin term of the GT operator also comes from the axial vector current and 
thus is closely related to the spin operators in Eq. i|l(j|) . Its is well established that most nuclear model 
calculations of GT strength require a reduction of the order of 20% in gA |3|j . Whether this quenching of 
gA should also be applied to a\ (the isovector neutralino-nucleon coupling constant) is unknown. Since 
there is no real guidance and magnetic moments obtained by Ressell and Dean agree well with experiment, 
it is very doubtful that any extra quenching of the spin matrix elements (or equivalently the coupling 
constants ao and ai) is desirable for these nuclei in the calculation of neutralino-nucleus scattering rates. 
Nonetheless, it is useful to keep these potential uncertainties in mind when calculating scattering rates 



Tables ITI lXIIII contain the fullest possible list of calculations of the nuclear zero- momentum spin prop- 
erties considered in the literature for detection of spin-coupled WIMPs. The tables are obtained on the 
basis of relevant tables for (S£) and (S£) from 0, S 0, S, 0, E3, EH E3, E3, EH EH IS3, H3| - The OGM 

results of Ellis and Flores given in |28l ISiA l4(l | in the form of A 2 J( J + 1) were recalculated into (Si) and 
checked. 



III. SUMMARY AND CONCLUSIONS 

There is continuous theoretical and experimental interest in existence of dark matter of the Universe. 
The best motivated non-baryonic dark matter candidates is the neutralino, the lightest supersymmetric 
particle. The motivation for supersymmetry arises naturally in modern theories of particle physics. In 
this work we discussed the spin-dependent interaction of neutralinos with odd mass nuclei. The nuclear 
structure plays an important role in determining the strength of the neutralino-nucleus cross section 
for this type of interaction. In the limit of zero momentum transfer the relevant physical quantities 
are the proton and neutron spin averages (S p ( n )), which have to be evaluated within a proper nuclear 
model. These values determine the event rate expected in a direct dark matter search experiment due 
to spin-dependent neutralino-nucleus interaction. In this work the calculation of spin-dependent matrix 
elements is reviewed. To our knowledge, a complete list of calculated spin matrix elements is presented 
for nuclei throughout the periodic table. We recall that only nuclei with an odd number of either protons 
or neutrons can have non-zero spin. 

As is manifested in this review, practically every known nuclear model has been employed for evaluation 



of the spin matrix elements. The results show that spin matrix elements depend in general rather 
sensitively on the details of the nuclear structure. The phenomenological ISPSM (independent single- 
particle shell model) is fairly accurate only for light nuclei with near-closed shells, but in general it 
tends to overestimate the spin matrix elements and is inadequate. It was confirmed by the calculations 
within the OGM (odd group model) and the EOGM (extended OGM) utilizing magnetic moments and 
mirror (3 decays. However, detailed shell model calculations found these phenomenological models to be 
inadequate. The odd group model and shell model treatments yielded good agreement for light nuclei, but 
as the atomic mass increases, there arouses a significant amount of configuration mixing not considered 
in the OGM. Unfortunately, the shell model calculations are difficult for most medium-heavy and heavy 
isotopes because of the size of the matrices involved. The situation is improving due to advances in 
computer power and storage. There is a hope to construct model spaces that contain most of the nuclear 
configurations that are likely to dominate the spin response of nuclei. For open-shell medium-heavy and 
heavy nuclei the methods of choice for calculation of spin matrix elements are the Interacting Boson- 
Fermion Model, the Theory of Finite Fermion Systems and the Quasi Tamm-Dancoff Approximation. 
The most reliable values of (S p ( n )) are considered to be those of the approach which reproduces well 
the experimental value of the magnetic moment for a given isotope. The magnetic moment is extremely 
important, as it is the observable most closely related to the neutralino- nucleus scattering matrix element 
and has traditionally been used as a benchmark for the calculation accuracy. 

There is an additional complication arising from the fact that the neutralino appears to be quite 
massive, perhaps heavier than 100 GeV. For such a heavy light supersymmetric particle and sufficiently 
heavy nuclei, the dependence of the nuclear matrix elements on the momentum transfer cannot be ignored. 
This affects the spin matrix elements. The calculations of the structure functions in the finite momentum 
approximation and the level of accuracy of these calculations are beyond the scope of this review. 

This work was supported in part by the VEGA Grant Agency of the Slovac Republic under contract 
No. 1/0249/03 and by the Russian Foundation for Basic Research (grant 02-02-04009). 



TABLE I: Zero momentum spin structure of light nuclei (A < 13) in different models. The measured magnetic 
moments used as input are enclosed in parentheses. 



H (Lj — 5*1/2) 


<s P ) 


(Sn> 


fx (in /in) 


ISPSM, Ellis-Flores [28, 39] 
OGM, Ellis-Flores [28, J59] 


1/2 
0.5 






2.793 

yZ. /yojcxp 


3 He (Lj = Si /a ) 


<s P ) 


(Sn> 


fx (in /i N ) 



ISPSM, Ellis-Flores [28., .22J 1/2 -1.913 

OGM, Engel-Vogel [29j 0.56 (-2.128) exp 

EOGM {g A /gv = 1), Engel-Vogel [29] -0.081 0.552 (-2.128) exp 

EOGM (g A /gv = 1-25), Engel-Vogel [2J -0.021 0.462 (-2.128) e x P 



7 Li (Lj = P3/2) 


<s P ) 


(Sn> 


fx (in fx N ) 


ISPSM, Ellis-Flores [28., .39] 
OGM, Engel-Vogel [29] 
SM, Pacheco-Strottman [14] 


1/2 

0.38 

0.497 






0.004 


3.793 
(3.256) cxp 


9 Be (Lj = P3/2) 


<s P ) 


(Sn> 


H (in /i N ) 


ISPSM, Ellis-Flores [28, 39] 
OGM, Engel-Vogel [29] 
SM, Pacheco-Strottman [14] 






0.007 


1/2 

0.31 

0.415 


-1.913 

(-1.178) oxp 


n B (Lj = P 3/2 ) 


<s P ) 


<s„> 


fx (in /i N ) 



ISPSM, Ellis-Flores [28, 39] 

OGM, Engel-Vogel [29j 

EOGM (g A /gv = 1), Engel-Vogel [29] 

EOGM (g A /gv = 1.25), Engel-Vogel [29] 

SM, Pacheco-Strottman [141 



1/2 





3.793 


0.264 





(2.689) cxp 


0.292 


0.006 


(2.689) oxp 


0.264 


0.034 


(2.689) cxp 


0.292 


0.008 





TABLE II: Zero momentum spin structure of light nuclei (11 < A < 21) in different models. 



13 C (Lj = P 1/2 ) 


<s P ) 


(Sn> 


(j, (in /i N ) 


ISPSM 





-0.167 


0.638 


OGM 





-0.183 


(0.702) cxp 


EOGM {g A /gv = 1), Engel-Vogel [29] 


-0.009 


-0.172 


(0.702) cxp 


EOGM (g A /gv = 1.25), Engel-Vogel [29] 


-0.026 


-0.155 


(0.702) cxp 


15 N (Lj = P 1/2 ) 


(Sp) 


(Sn> 


/U (in /ijv) 


ISPSM, Engel-Vogel [29] 


-0.167 





-0.264 


OGM, Engel-Vogel [22] 


-0.167 





(-0.283) cxp 


EOGM (g A /gv = 1), Engel-Vogel [29] 


-0.145 


0.037 


(-0.283) cxp 


EOGM (g A /gv = 1-25), Engel-Vogel [29] 


-0.127 


0.019 


(-0.283) cxp 


17 (Lj = D 5/2 ) 


<s P ) 


(Sn> 


IX (in /in) 


ISPSM, Ellis-Flores [28, 39] 





1/2 


-1.913 


OGM, Engel-Vogel [29J 





0.49 


(-1.894)exp 


EOGM (gA/ffv = 1), Engel-Vogel [29] 


-0.036 


0.508 


(-1.894) oxp 


EOGM (Wffv = 1.25), Engel-Vogel [29J 


0.019 


0.453 


(-1.894) cxp 


SM, Pacheco-Strottman [14] 





0.5 




19 F (Lj = S 1/2 ) 


<s P ) 


(Sn> 


fj, (in /xjv) 


ISPSM, Ellis-Flores [28, 39] 


1/2 





2.793 


OGM, Engel-Vogel [29] 


0.46 





(2.629) oxp 


EOGM {g A /gv = 1), Engel-Vogel [22] 


0.415 


-0.047 


(2.629) cxp 


EOGM (g A /gv = 1.25), Engel-Vogel [2£] 


0.368 


-0.001 


(2.629) cxp 


SM, Pacheco-Strottman [14] 


0.441 


-0.109 




SM, Divari et al. [5] 


0.4751 


-0.0087 


2.91 



TABLE III: Zero momentum spin structure of light nuclei (19 < A < 29) in different models. 



L Ne (Lj 



3/2J 



(Sp) 



(Sn> 



(j. (in /xjv) 



ISPSM 1/2 

OGM 0.173 

EOGM (g A /gv = 1), Engel-Vogel [29] 0.020 0.294 

EOGM (g A /gv = 1.25), Engel-Vogel [29] 0.047 0.2646 



5 Na (Lj 



■a/2) 



<S P ) 



(S„> 



-1.913 
(-0.662) cxp 
(-0.662) cxp 
(-0.662) cxp 



H (in fi N ) 



ISPSM 

SM, Ressell-Dean [6] 
OGM, Ressell-Dean [6] 
SM, Divari ar al. [5J 



1/2 
0.2477 
0.1566 
0.2477 




0.0198 

0.0 
0.0199 



5 Mg (Lj 



D, 



/2) 



<s P > 



(S.) 



3.793 
2.2196 

(2.218) cxp 
2.22 



(X (in fi N ) 



ISPSM 1/2 

OGM 0.223 

EOGM (g A /gv = 1), Engel-Vogel [29] 0.040 0.376 

EOGM (g A /gv = 1.25), Engel-Vogel [29] 0.073 0.343 



7 A1 (Lj = D 5/2 ) 



<S P ) 



(Sn> 



-1.913 

(-0.855) cxp 
(-0.855) cxp 
(-0.855) cxp 



H (in fi N ) 



ISPSM, Ellis-Flores [28., ,22] 1/2 4.793 

OGM, Engel-Vogel [29] 0.25 (3.642) cxp 

EOGM (g A /gv = 1), Engel-Vogel [22] 0.333 0.043 (3.642) cxp 

EOGM (g A /gv = 1.25), Engel-Vogel [2£] 0.304 0.072 (3.642) cxp 

SM, Engel et al. ]7j 0.3430 0.0296 3.584 



TABLE IV: Zero momentum spin structure of light nuclei (29 < A < 41) in different models. 

29 Si (Lj = S 1/2 ) (S p ) (S„) fi (in fi N ) 

ISPSM, Ellis-Flores [28, 39] 1/2 -1.913 

OGM, Engel-Vogel [29j 0.15 (-0.555) exp 

EOGM {g A /gv = 1), Engel-Vogel [29] 0.054 0.204 (-0.555) exp 

EOGM (g A /gv = 1.25), Engel-Vogel [29] 0.069 0.189 (-0.555) exp 

SM, Ressell et al. [g] -0.002 0.13 -0.50 

SM, Divari et al. [5] -0.0019 0.1334 -0.50 

31 P (Lj = 5 1/2 ) (S p ) (S„) M (m Mjv) 

ISPSM 0.5 2.793 

OGM 0.138 (1.132) exp 

EOGM (g A /gv = 1), Engel-Vogel [29] 0.181 0.032 (1.132) cxp 

EOGM {g A /gv = 1.25), Engel-Vogel [29] 0.166 0.047 (1.132) cxp 

35 C1 (Lj = D 3/2 ) (S p ) (SQ m (in /ijy) 

ISPSM, Ellis-Flores [28, 39] -0.3 0.13 

OGM, Engel-Vogel [22] -0.15 (0.822) cxp 

EOGM, Engel-Vogel [2£] -0.094 0.014 (0.822) cxp 

SM, Pacheco-Strottman [14] -0.059 -0.011 

SM, Ressell et al. [8] -0.051 -0.0088 

39 K [Lj = D 3/2 ) (S p ) (S„) m (in fi N ) 

ISPSM -0.3 0.324 

OGM -0.242 (0.391) cxp 

EOGM (g A /gv = 1.0) Engel-Vogel [29] -0.196 0.055 (0.391) cxp 

EOGM (g A /gv = 1-25) Engel-Vogel [29J -0.171 0.030 (0.391) cxp 

PT with Force I, Engel et al. [JJ -0.197 0.051 0.420 

PT with Force II, Engel et al. [7] -0.184 0.054 0.181 



TABLE V: Zero momentum spin structure of some nuclei (45 < A < 73) in different models. 



47 Ti (Lj = F 5/2 ) 


<s P > 


(Sn) 


/j, (in /ijv) 


ISPSM, Ellis-Flores [28, 39] 
OGM, Engel-Vogel [29j 






-0.357 
0.21 


1.367 
(-0.788) cxp 


49 Ti (Lj = F 7/2 ) 


<Sp> 


(Sn) 


;U (in /ijv) 


ISPSM, Ellis-Flores [2£, .3J 
OGM, Engel-Vogel [29] 






1/2 

0.29 


-1.913 

(-1.104) cxp 


51 V (Lj = F 7/2 ) 


<Sp> 


(Sn) 


/U (in fj, N ) 


ISPSM, Ellis-Flores [28, 39] 
OGM, Engel-Vogel [29] 


1/2 
0.36 






5.79 

(5.149) cxp 


55 Mn (Lj = F 5/2 ) 


(Sp> 


(Sn) 


/U (in /ijv) 


ISPSM, Ellis-Flores [£2, 40J 
OGM, Ellis-Flores [39, 40J 


-0.357 
0.264 






5.79 
(3.453) cxp 


59 Co (Lj = F 7/2 ) 


(Sp> 


(Sn) 


/U (in hn) 


ISPSM, Ellis-Flores [39, 40] 
OGM, Ellis-Flores [3£, 40J 


1/2 
0.25 






5.79 

(4.627) cxp 


67 Zn (Lj = F 5/2 ) 


<Sp> 


<s») 


^ (in /xjv) 


ISPSM, Ellis-Flores [28, 39j 
OGM, Engel-Vogel [29] 






-0.357 
-0.23 


1.367 

(0.875) cxp 


69 Ga (Lj = P 3/2 ) 


<Sp> 


(Sn) 


^ (in /xjv) 


ISPSM, Ellis-Flores [28, 39] 
OGM, Engel-Vogel [29] 


0.5 
0.11 






3.793 

(2.017) cxp 


71 Ga (Lj = P3/2) 


(Sp> 


(Sn) 


/i (in /ijv) 


ISPSM, Ellis-Flores [28, 39] 
OGM, Engel-Vogel [29J 


0.5 
0.23 






3.793 
(2.562) cxp 



TABLE VI: Zero momentum spin structure of some nuclei (71 < A < 95) in different models. 



J Ge (Lj 



G 



9/2 J 



(s P ) 



(Sn> 



fj, (in fi N ) 



ISPSM, Ellis-Flores [28, J39] 

OGM, Engel-Vogel [29J 

IBFM, Iachello et al. [15] and [8] 

IBFM (quenched), Iachello et al. [15] and [8] 

TFFS, Nikolaev-Klapdor-Kleingrothaus, [^OJ 

SM (small), Ressell et al. [8J 

SM (large), Ressell et al. [8] 

SM (large, quenched), Ressell et al. [8J 

"Hybrid" SM, Dimitrov et al. [9J 



ISPSM, Ellis-Flores [28,^9] 
OGM, Engel-Vogel [29] 
SM (large), Engel et al. [lfj] 
SM (small). Engel et al. [101 






0.5 


-1.913 





0.23 


(-0.879) cxp 


0.009 


0.469 


-1.785 


0.005 


0.245 


(-0.879) cxp 





0.34 


— 


0.005 


0.496 


-1.468 


0.011 


0.468 


-1.239 


0.009 


0.372 


(-0.879)exp 


0.030 


0.378 


-0.920 



75 As (Lj = P 3/2 ) 


<s P > 


(Sn) 


/i (in fi N ) 


ISPSM, Ellis-Flores [28,^9] 
OGM, Engel-Vogel [29] 


0.5 

-0.01 






3.793 
(1.439) cxp 


79 Br (Lj = P3/2) 


(S P > 


(Sn) 


fi (in fi N ) 


ISPSM, Ellis-Flores [2^,^9J 
OGM, Engel-Vogel [29J 


0.5 
0.13 






3.793 
(2.106) cxp 


81 Br (Lj = P3/2) 


<Sp> 


(Sn) 


fi (in fi N ) 


ISPSM, Ellis-Flores [28,^9] 
OGM, Engel-Vogel [29] 


0.5 

0.17 






3.793 

(2.271) cxp 


91 Zr (Lj = Z3 5/2 ) 


<s P > 


(Sn) 


IX (in hn) 


ISPSM, Ellis-Flores [2^,^9J 
OGM, Engel-Vogel [29J 






0.5 
0.34 


-1.913 

(-1.304) oxp 


93 Nb (Lj = G m ) 


<s P > 


(Sn) 


fx (in fi N ) 



0.5 





6.793 


0.36 





(6.171) cxp 


0.48 


0.04 


6.36 


0.46 


0.08 


5.88 



TABLE VII: Zero momentum spin structure of some nuclei (95 < A < 115) in different models. 



"Ru (Lj = D m ) 


<s P > 


(Sn> 


fi (in fi N ) 


ISPSM, Ellis-Flores [39, 40] 
OGM, Ellis-Flores [3^. 401 






1/2 
0.17 


-1.913 
(-0.6381) oxp 


101 Ru (Lj = D 5/2 ) 


<s P > 


(Sn> 


fi (in fj, N ) 


ISPSM, Ellis-Flores [39, 40] 
OGM, Ellis-Flores [39, 40] 






1/2 
0.19 


-1.913 
(-0.719) cxp 


107 Ag (Lj = P 1/2 ) 


<Sp> 


(S»> 


/U (in /ijv) 


ISPSM, Ellis-Flores [28, 39] 
OGM, Engel-Vogel [29J 


-0.167 
-0.13 






-0.264[-0.07] 
(-0.114) cxp 


109 Ag (Lj = P 1/a ) 


<Sp> 


(S»> 


jjb (in /ijv) 


ISPSM, Ellis-Flores [28, 39] 
OGM, Ellis-Flores [39, 40] 


-0.167 
-0.14 






-0.264[-0.07] 

(-0.131)exp 


m Cd (Lj = 5i /a ) 


<Sp> 


(S»> 


/z (in /ijv) 


ISPSM, Ellis-Flores [28, 39] 
OGM, Engel-Vogel [2£] 






1/2 
0.16 


-1.913 

(-0.595)exp 


113 Cd (Lj = S 1/2 ) 


<Sp> 


<s„) 


/U (in ^ijv) 



ISPSM, Ellis-Flores [28, 39] 1/2 

OGM, Engel-Vogel [29] 0.16 

IBFM, Iachello et al. [15J -0.001 0.488 

IBFM (quenched), Iachello et al. [15J -0.0 0.162 

TFFS, Nikolaev-Klapdor-Kleingrothaus, [30] 0.175 



5 Cd (Lj = Sua 



<Sp> 



(S»> 



-1.913 
(-0.622) cxp 



(i (in /njy) 



ISPSM, Nikolaev-Klapdor-Kleingrothaus, [30] 1/2 

IBFM, Iachello et al. [15J -0.001 0.488 

IBFM (quenched), Iachello et al. [15] -0.0 0.168 

TFFS, Nikolaev-Klapdor-Kleingrothaus, [30] 0.195 

OGM 0.169 



(-0.6388) e 



TABLE VIII: Zero momentum spin structure of heavy nuclei (114 < A < 125) in different models. Calculations 
of the magnetic moment using effective g-factors are given in curly brackets. 



115 Sn (Lj = S 1/2 ) 


(Sp) 


(Sn) 


(j, (in fi N ) 


ISPSM, Ellis-Flores [28., JJ9J 
OGM, Engel-Vogel [29] 






1/2 
0.24 


-1.913 
(-0.919)«xp 


117 Sn (Lj = Si /a ) 


<s P > 


<Sn> 


fj. (in fi N ) 


ISPSM, Ellis-Flores [28, 391 
OGM, Engel-Vogel [29J 






1/2 
0.126 


-1.913 

(-l.OOl)exp 


121 Sb (Lj = D 5/2 ) 


(Sp) 


(Sn) 


/i (in fi N ) 


ISPSM, Ellis-Flores [28., JJ9J 
OGM, Ellis-Flores [39, 40] 


0.5 

0.188 






4.793 
(3.363) oxp 


123 Sb (Lj = G 7/2 ) 


<s„) 


(Sn) 


/i (in fi N ) 


ISPSM, Ellis-Flores [28, J59] 
OGM, Ellis-Flores [g9_. 401 


-0.389 
-0.207 






1.717 
(2.550) cxp 


123 Te (Lj = S 1/2 ) 


(S„> 


(Sn) 


/i (in fi N ) 



ISPSM, Ellis-Flores [28,^9] 1/2 

IBFM, Iachello et al. [15] -0.000 0.491 

IBFM (quenched), Iachello et al. [15J -0.000 0.192 

TFFS, Nikolaev-Klapdor-Kleingrothaus, [^0J 0.21 

OGM 0.192 



5 Te (Lj = S 1/2 



(S P ) 



(Sn) 



ISPSM, Ellis-Flores [28,jB9] 1/2 

OGM, Engel-Vogel [2£] 0.0 0.23 

IBFM, Iachello et al. [15] -0.0008 0.499 

IBFM (quenched), Iachello et al. [15] -0.0004 0.231 

TFFS, Nikolaev-Klapdor-Kleingrothaus, [£0J 0.22 

SM (Bonn A), Ressell-Dean [6J 0.001 0.287 

SM (Nijmegen II), Ressell-Dean [6] -0.0003 0.323 



-1.913 



(-0.737). 



(j, (m fi N ) 



-1.913 

(-0.889) cxp 
(-0.889) cxp 
(-0.889) cxp 

-1.015 {-0.749} eff 
-1.134 {-0.824} off 



TABLE IX: Zero momentum spin structure of heavy nuclei (125 < A < 127) in different models. Calculations of 
the magnetic moment using effective <?-factors are given in curly brackets. 



T (Lj = D 5/2 



(Sp) 



(Sn) 



(j, (in (j, N ) 



ISPSM, Engel-Vogel [29] 1/2 

IBFM, Iachello et al. [15J 0.460 0.005 

IBFM (quenched), Iachello et al. [15] 0.159 0.002 

TFFS, Nikolaev-Klapdor-Kleingrothaus, [30] 0.18 

OGM 0.07 



7 l(Lj 



D 



5/2 J 



(Sp) 



(Sn) 



4.793 



(2.821) 



fj, (m )i N ) 



ISPSM, Ellis-Flores [^9., 40J 1/2 

OGM, Engel-Vogel [29] 0.07 

IBFM, Iachello et al. [15J 0.464 0.010 

IBFM (quenched), Iachello et al. [15] 0.154 0.003 

TFFS, Nikolaev-Klapdor-Kleingrothaus, [30] 0.15 

SM (Bonn A), Ressell-Dean [6] 0.309 0.075 

SM (Nijmegen II), Ressell-Dean [6J 0.354 0.064 



4.793 
(2.813) cxp 
(2.813) cxp 
(2.813) cxp 

2.775 {2.470} cff 
3.150 {2.7930} cff 



TABLE X: Zero momentum spin structure of heavy nuclei (128 < A < 133) in different models. Calculations of 
the magnetic moment using effective <?-factors are given in curly brackets. 

~ 29 Xe (Lj = S 1/2 ) (S„) (S w ) (j, (in fi N ) 

ISPSM, Ellis-Flores [28, 39] 1/2 -1.913 

OGM, Engel-Vogel [29] 0.0 0.2 (-0.778) OX p 

IBFM, Iachello et al. [15J -0.000 0.430 (-0.778) oxp 

1BFM (quenched), Iachello et al. [15J -0.000 0.200 (-0.788) exp 

TFFS, Nikolaev-Klapdor-Kleingrothaus, [30] 0.25 

SM (Bonn A), Ressell-Dean [6] 0.028 0.359 -0.983 {-0.634} eff 



SM (Nijmegen II), Ressell-Dean [6] 0.0128 0.300 -0.701{-0.379} 



eff 



131 Xe (Lj = £> 3/2 ) (S p ) (gn) (M (in fi N ) 

ISPSM, Ellis-Flores [28, 39] -0.3 1.148 

OGM, Engel-Vogel [29] 0.0 -0.18 (0.692) cxp 

IBFM, Iachello et al. [1.5] 0.000 -0.280 (0.692) exp 

IBFM (quenched), Iachello et al. [15J 0.000 -0.168 (0.692) exp 

TFFS, Nikolaev-Klapdor-Kleingrothaus, [30] -0.186 

SM (Bonn A), Ressell-Dean [6] -0.009 -0.227 0.980 {0.637} cff 

SM (Nijmegen II), Ressell-Dean [6J -0.012 -0.217 0.979 {0.347} c ff 

QTDA, Engel [12] -0.041 -0.236 0.70 

133 Xe (Lj = D 3/2 ) (S p ) (S n ) ii (in /i N ) 

ISPSM, Ellis-Flores [28, 39] -0.3 1.148 

IBFM, Iachello et al. [15J 0.000 -0.257 

IBFM (quenched), Iachello et al. [15J 0.000 -0.176 

TFFS, Nikolaev-Klapdor-Kleingrothaus, [30] -0.201 

OGM 0.0 -0.213 (0.813) exp 



TABLE XI: Zero momentum spin structure of heavy nuclei (133 < A < 141) in different models. 



133 Cs (Lj = G r /a) 






<s P > 


(Sn) 


fi (in ^ijv) 


ISPSM, Ellis-Flores [28, J39] 






-0.389 





1.717 


OGM, Engel-Vogel [29J 






-0.20 





(2.582) cxp 


IBFM, Iachello et al. [15] 






-0.370 


0.003 




IBFM (quenched), Iachello et al. 


[15] 




-0.225 


0.002 




TFFS, Nikolaev-Klapdor-Kleing: 


rothaus, 


m 


-0.230 






135 Cs (L., = G r /a) 






<s P > 


(S») 


H (in ^ijv) 


ISPSM, Ellis-Flores [2^,^9J 






-0.389 





1.717 


OGM 






-0.167 





(2.734) cxp 


IBFM, Iachello et al. [15] 






-0.373 


0.002 




IBFM (quenched), Iachello et al. 


ri5] 




-0.201 


0.001 




TFFS, Nikolaev-Klapdor-Kleing: 


rothaus, 


[30] 


-0.199 






135 Ba (Lj = D 3/2 ) 






<s P > 


(Sn) 


fi (in fi N ) 


ISPSM, Ellis-Flores [28,^9] 









-0.30 


1.148 


OGM 









-0.219 


(0.838) cxp 


IBFM, Iachello et al. [15J 






-0.007 


-0.226 




IBFM (quenched), Iachello et al. 


[15] 




-0.004 


-0.145 




TFFS, Nikolaev-Klapdor-Kleingrothaus, 


[30] 




-0.18 




137 La (Lj = G 7/2 ) 






<s P > 


(Sn) 


fi (in ^ijv) 


ISPSM, Ellis-Flores [2^,^9J 






-0.389 





1.717 


OGM 






-0.176 





(2.695) cxp 


IBFM, Iachello et al. [15] 






-0.386 


0.0006 




IBFM (quenched), Iachello et al. 


[15] 




-0.212 


0.0003 




139 La (Lj = G 7/2 ) 






<s P > 


(Sn) 


fi (in ^ijv) 


ISPSM, Ellis-Flores [2^,^9J 






-0.389 





1.717 


OGM, Engel-Vogel [29] 






-0.16 





(2.783) cxp 



TABLE XII: Zero momentum spin structure of heavy nuclei (143 < A < 205) in different models. 



155 Gd (Lj = P 3/a ) 


<s P > 


(Sn) 


H (in fi N ) 


ISPSM, Ellis-Flores [39, 40] 
OGM, Ellis-Flores [39, 40] 






0.5 
0.07 


-1.913 

(-0.259)exp 


157 Gd (Lj = P 3/2 ) 


<s P > 


(Sn) 


H (in fi N ) 


ISPSM, Ellis-Flores [39, 40] 
OGM, Ellis-Flores [39, 401 






0.5 
0.09 


-1.913 
(-0.340) cxp 


183 W (Lj = P 1/2 ) 


<s P > 


(Sn> 


/U (in /ijv) 


ISPSM, Ellis-Flores [39, 40] 
OGM, Ellis-Flores [39.401 






-0.17 
-0.03 


0.638 

(0.118) cxp 


191 Ir (Lj = D 3/2 ) 


<s P > 


(Sn> 


jjb (in /ijv) 


ISPSM, Ellis-Flores [39, 40] 
OGM, Ellis-Flores [39, 40] 


-0.30 
-0.295 






0.148 

(0.151)exp 


193 Ir (Lj = D 3/2 ) 


<s P > 


(Sn> 


/z (in /ijv) 


ISPSM, Ellis-Flores [39, 40] 
OGM, Ellis-Flores [£9., 40J 


-0.30 
-0.292 






0.148 

(0.164)exp 


199 Hg (Lj = P 1/2 ) 


<s P > 


<s„) 


fi (in fj, N ) 


ISPSM, Ellis-Flores [28, 39] 
OGM, Engel-Vogel [29J 






-0.17 
-0.13 


0.638 

(0.506) cxp 


201 Hg (Lj = P 3/2 ) 


<s P > 


<s„) 


H (in fi N ) 


ISPSM, Ellis-Flores [39, 40] 
OGM, Ellis-Flores [39, 401 






0.5 
0.146 


-1.913 
(-0.560) cxp 


203 T1 (Lj = Si /a ) 


<s P > 


(Sn> 


/U (in /ijv) 


ISPSM, Ellis-Flores [28,, 39] 
OGM, Engel-Vogel [29] 


0.50 
0.24 






2.793 

(1.662)exp 


205 T1 (Lj = Si/ a ) 


<s P > 


(S»> 


/z (in /ijv) 


ISPSM, Ellis-Flores [28, 39] 
OGM, Engel-Vogel [29] 


0.50 
0.25 






2.793 

(1.638)exp 



TABLE XIII: Zero momentum spin structure of heavy nuclei (A < 209) in different models. 



207 Pb (Lj = P 1/2 ) 


<s P ) 


(Sn> 


/U (in /ijv) 


ISPSM, Ellis-Flores [28, 39J 

OGM, Engel-Vogel [29J 

SM, Kosmas-Vergados [3, 13] 






-0.010 


-0.17 

-0.15 

-0.149 


0.638 

(0.593) cxp 


Bi (Lj — H9/2) 


(Sp) 


(Sn> 


^ (in /ijv) 


ISPSM, Ellis-Flores [^9., 40J 
OGM, Ellis-Flores [39, 40] 


-0.41 
-0.085 






2.63 

(4.111)exp 



IV. APPENDIX A 

A. Elements of nuclear structure calculations 

The transverse electric T el5 (q) and longitudinal C 5 (q) multipole projections of the axial vector current 
operator as well as the scalar function Cl (q) are given by la, \12L l24j : 

r? ( a ) - 1 W^£ + Q l m ' T 3 \ 

[Vi + lM LiI , + i(gr-) + vTAfi.i-iCgfi)] , 

C L (q) = ]T CfMmWLifi), c (q)=J2 c oJo(m)Yo(h), (22) 

i, nucleons i 

where ML,L'(QTi) — j l' (q r i)[YL' {vi)<?i\ L . In the limit of zero momentum transfer S^ D (q) reduces to 

i 

= -^|(ao + ai)(7V||S p ||7V) + ( ao - ai )(iV||S„||JV)| 2 (23) 

47T 

1(2J+1)( J + !),__ /ArlQ |An , _ /Ar)Q 2 



7T J 

2J+1 



|ap<JV|S p |JV)+a n (iV|S n |JV>|' 



-J(J +1)A 2 , (24) 

7T 
Wlm_fi_ J-/J._S_IAA n._/SA 

with A = 



7T 

(A^|a p Sp + a n S n |./V) _ a p (S. p ) a n (S n ) 



In accordance with convention the Z components of the angular momentum and spin operators are 
evaluated in the maximal Mj state, e.g. (S) = (JV|S|iV) = (J,Mj = J\S z \J,Mj = J). 

In the ISPSM only the last odd nucleon contributes to the spin and the angular momentum of the 
nucleus. In this limit 

(Sft = W' + "-;ff' + " + « , (25, 

where J a and La are the single-particle total and angular momenta. They are deduced from the measured 
nuclear angular momentum and the parity. 

B. Nucleon spin structure 

To evaluate the spin content of the nucleon one needs the matrix element of the effective quark axial- 
vector current J M = qj^jsq in the nucleon |l6(. These matrix elements 

((P,n)|fvy*g|(p,n)) = 2»<f»">Ag<»"> (26) 

are proportional to the spin of the neutron (proton or neutron), s{? . The quantities Aq^' 71 ' are usually 
extracted from the data obtained in polarized lepton-nucleon deep inelastic scattering. Uncertainties 
in the experimentally determined values for the quantities Aq can lead to significant variations in the 
WIMP-nucleon axial- vector coupling, and therefore to the predicted rates for detection of WIMPs which 
have primarily spin couplings to nuclei |l6( . With definition (|26|) the effective spin-dependent interaction 
of neutralinos with the nucleon has the form 

£ spin = 2x7"75X ns^n ^ A q Aq^. (27) 

g— u,d,s 

Recent global QCD analysis for the g\ structure functions |4l|, including 0{p? a ) corrections, corresponds 
to the following values of spin nucleon parameters [42J 

A (P) _ A W = q 7g ± q Q2) A (P) = A (n) = _ Q 4g ± q 02) A (p) _ A („) = _ Q 15 ± Q 02 



X 



X 




FIG. 1: Spin-dependent elastic scattering of neutralinos from quarks. 

X XX X 

q 

q q ^ ^V q 



H,h 




FIG. 2: Spin-independent or scalar (tree level) elastic scattering of neutralinos from quarks. 



C. Effective neutralino-quark Lagrangian 

The axial-vector and scalar interaction of a neutralino with a quark q is given by 

£eff = A ' X7m75X • <Z7 M 75<7 + C g ■ XX ' QQ + 0(l/m~). 

The terms with vector and pseudoscalar quark currents are omitted being negligible in the case of non- 
relativistic DM neutralinos with typical velocities v x ~ 10~ 3 c. The Feynman diagrams which give rise to 
the effective neutralino-quark axial-vector couplings 



s%a 



9 



4M^ 



Af? A -Af?: 



'14 - ,V 13 T 

— l a 



Mw 



Ml 



m q2 - ( m X + m q) 

!!i P 2 ( 



r (sin 2 



si - ("be + m i) : 



->q 9qL + COS V q (p qR ) 



+ 



:(cos 2 8 q 



& 



4 q \ m-j — (m x + m q ) 2 m? 2 — (m x + m q ) 2 
^ M w P q sin20 g T 3 (A/i 2 - tan^Mi) 



sin 2 q <t> 2 qR ) 



m |l ~~ ( m X + m q) 2 m2 q 2 ~ ( m X + m q) 2 



are shown in Fig.^ The first term in A q comes from Z° exchange, and the other terms come from squark 
exchanges. The Feynman diagrams which give rise to the effective neutralino-quark scalar couplings 



C a 



m q g* 

M w 4 



Fi 



i-i 






H„ 



•H 



( m q p 2 _ M w 
q 



sin 29 n 



w 

sin 29 „ 



m„ 



4>qL 4>qR 



51 - ( m x + m q) 2 m q2 - ( m x + m q) 2 



where 



P„ 



cos 2 9 q 4> qL - sin 2 6 q <j> qR cos 2 9 q <j) qR - sin 2 9 q 4> qL 



m \\ - ( m X + m q) 2 



l h -(m x +m q ) 2 



F h = (A/12 — Ml tan 0vp) (A/14 cosuh + A/13 sina_fj), 
Fh — (A/12 —A/11 tan Ow) (A/14 sin an —A/13 cos uh), 



,1 ,cosa# 

^2 + ^lmT 



,1 ,sina ff 

y 2 ' cosp 



4>qL = M2T3 + A/ii (Q-T 3 ) tan %, giJ = tan% Q A/ii, 



A = 



•2 3; sin/3 



-1 Tl AA 13 
•2 3; cos/3' 



are shown in Fig. EI The importance of these scalar spin-independent contribution was found by K.Griest 
in EH. 



D. SUSY particle spectrum 

For completeness, we collect here formulas for masses of the SUSY particles in the MSSM. There are 
four Higgs bosons — neutral CP-odd (A), CP-even (H,h), charged (H )■ The CP-even Higgs boson 
mass matrix has the form: 



Pn P12 
Pi 2 P22 



1 / tan/3 -1 

2 \ -1 cot/3 

1 /cot/3 -1 

2 I — 1 tan/3 



M% sin 2/3 



m z sin 2/3 + u> 



An A 
A12 A 2 



12 



#n = 



sin 2/3 ?n 



(^ + Mi tan/3) -fa, A u , 



P. 



sin 2/3 / , 



22 



(m^ tan/3 



a -)+ojA 22 , 



tan/3' 



#12 - #21 - - 



sin 2/3 . , , ,, . 



For example, An which includes loop corrections is 



A n 



<n m ii m b2 , 2A 6 (A-/itan/3) ] "^ 



-(In 



°/3 

m-j . Aj (At — /i tan /3) 



2 2 

m- — m~ 

61 fc 2 



In- 



HIT 



-) 



/--ov"o f- — "f/\2 ,-2 2 \ . "H 1 

( — Z3 —2 ) 9{m~ b ,m~ b J + — ( 






tan/3' -.2 



%3 



'/J 



rri 1 +m 2 



51 *2 



^ = 1K 2 2 ci — cos 8, si = sin 8, q(mi.mi) — 2 \ Jin— 7. The diagonalization of the above 

matrix gives the Higgs boson masses mji,h- 

A,h = \{Hii + P22 ± yJ(H n + P 22 ) 2 - 4(P n P 22 - P^)}, 



HIT 



m 2 H± 



l w 



M A + wA ch . 



Here ra/f± is the charged Higgs boson mass in the one-loop approximation. The mixing angle a H is 
obtained from 



sin 2a w = — 7 



2H l2 



cos2a H = 



Til Tjl 

-"11 - ^22 



"H° 



■H" 



The neutralino mass matrix in the basis (P, W , P", P 2 ) has the form 



M, 



I Mi — Mz cos 8 sinOw Mz sin /3 sin 9w \ 

M 2 M^ cos 8 cos #w — M z sin 8 cos 9w 

-Mz cos 8 sinOw Mz cos 8 cos 8\y — /z 

y M z sin /3 sin 9w —Mzs\n(3cos8w — A* / 



The diagonalization gives mass eigenstates (4 neutralinos) : 

Xl (m Xi ) = Mi P + M 2 H> 3 + M 3 P° + M 4 P 2 ° • 
The lightest (LSP) x — Xi is the best DM candidate. The chargino mass term is 



(w-,Sr) 



M 2 V2M w sia/3 \ ( W A 

V2M w cos(3 \x ) \ H% 



h.c. 



The diagonalization U*My±V 1< = diag(Af~± , M~± ) gives charged mass eigenstates 



U a W-+U a H-, x + = VaW + + V l2 H^ 



with masses 



M 2 



Xl,2 



[Ml + ^ + 2M 2 vT 



T V (M| - a* 2 ) 2 + 4M^ cos 2 2/3 + 4M 2 / (Af | + M 2 + 2M 2Ai sin 2/3) 



The sfermion mass matrices M~, M 2 and A4? have the form: 

t 7 b T 



M\ = 



Mi = 
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Ml 



m| + m 2 + i(4M 2 / - M|) cos 2/3 
m t (A t -^cot/3) 



m t (A t - jucot/?) 



2 + m 2 - |(M2, -M|) cos 2/3 



m 2 - \{2M^ + M§) cos 2/3 mfc(A 6 - //tan/3) 

m^ + m 2 + KM 2 , - Af |) cos 2/3 



TOfc(ylb — /itan/3) 



m\ + m 2 - ±( 2M w - M z) cos 2/3 
m T (A T — /ttan/3) 



m T (A T — /xtan/3) 
m\ + m 2 + (Af^ - Af|) cos 2/3 



It is worth noting that these masses as well as the above-mentioned couplings of neutralino-quark 
interactions A q and C q are functions of the common set of SUSY parameters like, for example, tan/3, 
Ma, /i, A q , etc. The set of parameters allows one to describe observables at the highest and lowerst 
energies coherently and simultaneously. 
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